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ABSTRACT 

We report the discovery of thermal X-ray emission from the youngest Galactic supernova remnant Gl. 9+0.3, 
from a 237-ks Chandra observation. We detect strong Ka lines of Si, S, Ar, Ca, and Fe. In addition, we 
detect a 4.1 keV line with 99.971% confidence which we attribute to 44 Sc, produced by electron capture from 
44 Ti. Combining the data with our earlier Chandra observation allows us to detect the line in two regions 
independently. For a remnant age of 100 yr, our measured total line strength indicates synthesis of (1 -7) x 
10~ 5 M Q of Ti, in the range predicted for both Type la and core-collapse supernovae, but somewhat smaller 
than the 2 x 10~ 4 M Q reported for Cas A. The line spectrum indicates supersolar abundances. The Fe emission 
has a width of about 28,000 km s -1 , consistent with an age of ~ 100 yr and with the inferred mean shock 
velocity of 14,000 km s -1 deduced assuming a distance of 8.5 kpc. Most thermal emission comes from regions 
of lower X-ray but higher radio surface brightness. Deeper observations should allow more detailed spatial 
mapping of 44 Sc, with significant implications for models of nucleosynthesis in Type la supernovae. 
Subject headings: ISM: individual objects (Gl.9+0.3) — ISM: supernova remnants — nuclear reactions, nu- 
cleosynthesis, abundances — X-rays: ISM 



1. INTRODUCTION 

The supe rnova remnant (SNR) Gl .9+0. 3 has an age of orde r 
100 years (Rey nolds et al.|2008, Pa"peTij |Green et al.|2008) . 
Its integrated X-ray spectrum is well described by a model 
of synchrotron emission from a power-law distribution with 
an exponential cutoff (XSPEC model srcut), with rolloff 
frequency hv TO \\ = 2.2 keV and a very high absorbin g col- 
umn density Nh- 5x 10 22 cm" 2 ( jReynolds et al. 2009 Paper 
II). The high column suggests a distance of order that to the 
Galactic Center; we adopt a nominal distance of 8.5 kpc, at 
which the mean expansion rate is v s ~ 14,000 km s" 1 (Paper 
I). No thermal emission is apparent in the integrated spectrum 
based on short (50 ks in duration) Chandra observations in 
early 2007. Without the detection of thermal emission, basic, 
crucial information such as the supernova (SN) type, distance, 
and elemental abundances cannot be obtained. Its detection 
was one of the goals of a much longer Chandra observation. 



2. OBSERVATIONS 

We re-observed Gl.9+0.3 for 237 ks with Chandra in four 
observations between July 13 and 26, 2009, using the ACIS- 
S CCD camera (S3 chip). We checked aspect correction 
and created new level-1 event files appropriate for VFAINT 
mode. No flares occurred during the observation. CTI 
correction was applied and calibration was performed using 
CALDB version 3.4.0. Finally, the datasets were merged and 
weighted response files created. We extracted spectra using 
the specextract script. We obtained about 40,000 source 
counts. 

Figure [T] shows th e 2009 image, smooth e d with the spatio 
spectral method of Krishnamurthy et al. (2010). Figure " 
compares the radio and X-ray morphologies. The marked an 
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Figure 1. Chandra image of Gl.9+0.3. Red, 1 - 3 keV; green, 3 - 4.5 keV; 
blue, 4.5 - 7.5 keV. Image size 127" x 121". 

surprising difference between them, particularly strong radio 
but weak X-rays in the north, suggests that the radio image 
may hold a clue to the location of thermal emission. We ex- 
tracted the spectrum from the radio-bright region shown on 
Figure [2j it is shown in Figure [3] Lines are apparent. We 
found lines as well in the central region shown in Figure |2j 
and analyzed the regions separately. 

3. SPECTRAL ANALYSIS 

We combined the 2007 and 2009 datasets for a total ex- 
posure of 286 ks. We modeled the background rather than 
subtracting it, and used Markov chain Monte Carlo (MCMC) 
methods as implemented in the PyMC software package (Patil 
et al. 2010| > to determine best p arameter values and error 
ranges (e.g., |van Dyk et a l. 2001). MCMC methods require 
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in the emission of X-ray and 7-ray lines in very you ng rem 
nants ( 44 Ti decays with a mean life of 85.0 ±0.4 yr, Ahmad 



et al. 2006| l. This decay commences via an electron capture 



Figure 2. Top: 1.4 GHz radio image (Green et al. 2010, in preparation). 
Resolution 2.3" X 1.4". Bottom: X-ray image showing regions from which 
spectra below were extracted. Blue: North rim. Red: center. 

the specification of priors on parameters to be determined (as 
described below). 

The spectrum of the northern, radio-bright rim is shown in 
the upper panel of Figure [3] Spectral lines typical of strongly 
underionized plasma are apparent (such plasma is expected in 
Gl.9+0.3 because of its youth and the low density of the am- 
bient ISM). SN 1006 has a very simil ar X-ray spectrum; in 
the same spectral range of Figure 0] Yamaguchi et al. (2008) 
find prominent Ka lines of abundant elements such as Si, S, 
Ar, Ca, and Fe, with line centroids at 1.815 keV, 2.36 keV, 
3.01 keV, 3.69 keV, and 6.43 keV. (As in SN 1006, O, Ne, and 
Mg lines might also be present at lower energies, but cannot 
be seen because of the high absorption.) In addition to these 
lines produced in hot shocked plasma, the radioactive decay 
of 44 Ti to 44 Sc and finally to the stable isotope 44 Ca will result 



to 44 Sc, leaving a K-shell vacancy followed rapidly either by 
Auger decay or by emission of a fluorescence photon of en- 
ergy 4.09 keV ( the yield is 0. 172 photons per each 44 Ti decay; 
Be et al. 2004). Nuclear de-excitation Sc 7-ray lines at 78.4 



and 67.9 keV are also emitted, followed by Ca gamma rays at 
1.157 MeV (mean life of 44 Sc is 5.4 hr). An inspection of the 
radio-bright rim spectrum (Figure [3]l reveals the presence of 
Si, S, Ar, and Fe Ka lines, and a broad feature near 4 keV that 
may be a blend of Doppler-broadened Ca and Sc Ka lines. 
Lines are generally weaker in the low-surface brightness inte- 
rior (lower panel of Figure |3J, with Ar and Sc lines being the 
most prominent. 

We modeled the spectra of the northern, radio-bright rim 
and the faint interior with an absorbed power law plus emis- 
sion lines of Si, S, Ar, Ca, Sc, and Fe. This simple model 
does not account for dust scattering and does not separate the 
underlying continuum into thermal and nonthermal compo- 
nents, but it suffices for the determination of line strengths, 
centroids, and widths. We used a normal (Gaussian) prior 
for the absorbing column density Nh, with mean (standard 
deviation) of 6.89(0.11) x 10 22 cm" 2 , based on our multire- 
gion spectral fit without dust scattering to the 2007 data (Pa- 
per II). (Solar abundances of the absorbing ISM are those of 
Grevesse & Sauval 1998; fits with dust scattering resulted in 
Nh lower by 25%.) Noninformative, uniform and logarithmic 
priors were assumed for the power-law index T and the (un- 
absorbed) 5-10 keV continuum flux .F5-10 kev> respectively. 

We used normal priors for line energies, setting mean ther- 



mal line ene rgies equal to the values measured by Yamaguchi 
|et~aLl ( |2008] > for SN 1006, and to 4.09 keV for the Sc line. 
In view of the 14,000 km s" 1 blast wave speed, significant 
bulk Doppler shifts are possible, so we chose a large (<r = 10 4 
km s" 1 ) width for these priors. (For numerical stability, these 
normal priors were truncated to include only a finite range 
in line energies; we verified that our results are not affected 
by this procedure.) Thermal lines in a young remnant arise 
in a fast-moving, shocked shell bounded by forward and re- 
verse shocks. An optically and geometrically thin shell ex- 
panding with velocity v s h e it produces flat-topped lines with 
Doppler widths of 2v S heii- We assumed flat-topped profiles 
with the same (but unknown) Doppler width for all thermal 
lines. (Therma| broadening is likely of only modest-impor- 
tance, (kTj/mj) 1 ~ 0.2 v S heii, based on models of Dwarkadas 
|& Chevalier| ([1998) with exponential ejecta density pronies. 
But it may still be appreciable in off-center locations such as 
the north rim, where bulk radial motions contribute less to 
the line broadening.) A truncated normal prior was assumed 
for v s h e ii, with mean of 14,000 km s" 1 and la width of 5000 
km s" 1 , extending from to 50,000 km s" 1 . The Sc line was 
modeled by a Gaussian with width cr~ c ; we assumed a half- 
normar] prior with tr = 0.15 keV for af e . These priors for 
line wraths exclude very large (> 50,000 km s" 1 ) widths, but 
otherwise provide weak constraints. 

Some constraints on the Sc line strength are provided by 
IBIS/ISGRI onboard INTEGRAL. The 68 an d 78 keV lines 
have been detected by IBI S/ISGRI in Cas A (|Renaud et al.| 
[2006l > but not in Gl.9+0.3 ( |Renaud et al.|2009| ). We use the 
Cas A detection to bound priors. We chose a weakly infor- 
mative gamma prior for the Sc line strength, defined as the 

5 normal distribution with mean limited to positive domain 
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expected line counts in a 286 ks exposure with Chandra, with 
shape parameter a = 1.2 and scale parameter (3 = 0.01 (see 
|van Dyk et al. (2001 ) for discussion of the use of gamma pri- 
ors in modeling emission lines in X-ray spectra). Its mean of 
a/ (3= 120 counts is comparable to the square root of the vari- 
ance, equal to (a//? 2 ) 1 ' 2 =110 counts. This prior disfavors 
a stro ng Sc line that was sea rched for but not found in Cas 
A (Theiling & Leising 2006 ), but otherwise provides rather 
weak constraints. We assumed the same weakly informative 
prior on line strength for thermal lines of Si, S, Ar, Ca, and 
Fe. 

The sky background was determined by a fit to the back- 
ground spectrum extracted from a large area on the S3 chip. 
The sky background was modeled by two absorbed power- 
laws, while the particle background model involved a combi- 
nation of power-laws with exponential cutoffs and narrow flu- 
orescent lines. We allowed for spatial variations in the particle 
background, with a logarithmic prior imposed on the particle 
normalization. 

Table [Tlcontains results of spectral fits; models are plotted 
in Figure [3] Except for fluxes, values quoted in Table [T] are 
standard means of the MCMC draws. We used a geometric 
mean for ^ e y, and mean line fluxes F w ere evaluated as 
((F 1 / 2 )) 2 (according to recommendations by van Dyk et al. 



2001 1. The 90% confidence intervals were computed using 



the 0.05 and 0.95 quantiles of the draws. The Sc line width is 
poorly constrained, so we provide only upper limits based on 
the 90% highest probability density interval for of c . 

The Sc line and thermal lines of all abundant elements have 
appreciable strengths in the spectrum of the northern rim. 
Thermal lines are very broad, with a full width of 28,000 km 
s" 1 , consistent with the estimated shock velocity of 14,000 km 
s -1 . The Fe Ka line provides the strongest constraints on line 
widths, but Doppler broadening is important for other lines 
as well. Photon statistics are insufficient to study variations 
in Doppler widths between different elements. This includes 
Sc, with a width that does not seem to be different from ther- 
mal line widths. But Sc forms a line blend with Ca, so errors 
are particularly large for these elements. In the interior, ther- 
mal lines are much weaker; an inspection of the confidence 
intervals in Table 1 suggests that only Ar may have been un- 
ambiguously detected. Notwithstanding this decline in the 
overall strength of thermal lines, the Sc line remains strong 
(it is only a factor of 2 weaker than in the northern rim). The 
prominence of Sc in the interior is consistent with an origin in 
the unshocked ejecta. 

In both spectral regions independently, we examined the 
significance of the S c line detectio n using a likelihood ratio 
test as described in Protassov et al.| ( |2002| >. The null model 
consisted of an absorbed power law plus Si, S, Ar, Ca and 
Fe lines described by flat line profiles with equal (but un- 
known) Doppler widths, together with the same priors that 
we used for our spectral fits. The alternative model is ob- 
tained by adding an additional line at 4.09 keV to the null 
model. Large samples from the posterior distribution of the 
null model were obtained using MCMC simulations. From 
these samples we simulated synthetic datasets, fit each dataset 
with the null model and the alternative model, and computed 
the likelihood ratio statistic. The posterior predictive p-values 
for the null model are 0.0098 and 0.030 for the northern rim 
and interior, respectively. The p-value is equal to 0.00029 for 
the null model to be valid in both of these regions. We can 
therefore reject the null model, i.e., claim the detection of the 
4.1 keV line, at 99.971% significance. 
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Figure 3. Top: Spectrum of N rim (radio-bright region) shown in Figure[2] 
Lines of Ca and Sc blend together because of the large Doppler widths. Bot- 
tom: Spectrum of interior region. In both cases, background has been mod- 
eled rather than subtracted; source and background models are described in 
the text. The two black lines are the background and source models; the red 
line is the total. 

The total Sc line flux was derived by combining the spectra 
shown in Figure[3j and fitting the combined dataset separately 
with the same spectral model as before. We obtain a total Sc 
line flux of 1.2 x 10" 6 ph cm s , with a 95% confidence in- 
terval of (0.35,2.4) x 10~ 6 ph cm -2 s" 1 . The column density of 
6.9 x 10 22 cm -2 implies that X-rays of 4. 1 keV are attenuated 
by about 33% due both to absorption and scattering. The Sc 
line fluxes have not been corrected for this, but the correction 
has been made for the 44 Ti masses that are discussed next. 

We expect 0.172 fluorescence photons per 44 Sc atom. Thus 
the line flux, along with the source age and distance, gives 
directly the amount of 44 Ti synthesized in the explosion. In- 
cluding the branching ratio for the emission of the 1 . 157 MeV 
photon, we expect the flux of X-ray photons to be related to 
that in 1.157 MeV photons by F x = 0.\5F 1 . The gamma-ray 
flux is related to the mass of 44 Ti by 



F 7 = 3.6x 10" 



D 



M(Ti) 

10" 5 M Q V8.5kpc 



-t/r 



ph cm 2 s 1 



with the mean life r = 85 yr. (If the Ti is ionized beyond 



the He-like state, the effective lifetime is longer (Mochizuki 
et al.|1999) , but the low ionization age of Fe implied by the 



line centroid reported in Table [TJmeans that this effect should 
be negligible.) We used these relations with an age of 100 yr 
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to obtain the following 44 Ti masses: 2.8(1.1,5.1) x 10" 5 M 
in the north, 1.3(0.36,2.6) x 10" 5 M Q in the interior, and 
3.3(0.95,6.5) x 10" 5 M Q for the total. (Errors are 90% con- 
fidence limits for the two subregions, but 95% limits for the 
total.) These masses have been corrected for absorption and 
scattering by being multiplied by a factor of 3/2. For an age 
of 140 yr (corresponding to no deceleration at all), the masses 
would be larger by an additional factor of 1 .54. 

We also modeled the spectrum of the northern, radio- 
bright rim with a plane shock non-equilibrium ioniza- 
tion (NEI) model (vp shock in XSPEC). We find N H = 
6.88(6.72,7.03) x 10 22 cm" 2 ; kT e = 3.6 ±0.4 keV; and ion- 
ization timescale r = n e t = 1.4(0.78,3.1) x 10 9 cm" 3 s. This 
extremely low ionization age is very unusual for SNRs, but 
appropriate for an unprecedentedly young object. Abun- 
dances have substantial errors, but solar abundances are ruled 
out: compared to solar, Si (set equal to S) = 3.4 (2.0, 4.6), 
Ar = 17 (2.8, 37), and Fe = 4.1(2.4, 5.7). A blueshift of 
line centers is required, v r = -4,300(-7,600,-l,300) km s" 1 . 
We find a FWHM of 26,500(17,600,34,600) km s" 1 . With 
n e ~ 2t/100 yr= 0.7 cm" 3 , the estimated shocked ejecta mass 
is only ~ O.O16M , with comparable masses of intermediate- 
mass (Si, S, and Ar) elements and Fe. This is a quarter of the 
total shocked ambient gas mass of 0.06M Q (for a preshock 
density of ~ 0.04 cm" 3 and a blast wave radius of 2.2 pc; Pa- 
per I). Impact of the ejecta with a denser than average ambient 
medium in the north may explain the relatively high shocked 
ejecta mass and density, and the strong radio emission there. 

4. DISCUSSION 
We summarize our results as follows: 

1 . The radio-bright N rim region shows strong Ka lines of 
Si, S, Ar, Ca, and Fe (much weaker thermal emission is 
also present in the shell interior). 

2. The spectral lines are both broadened and shifted. 
Widths of thermal lines are 28 , 000 ± 5 , 000 km s" 1 with 
blueshifts of about 4,000 km s" 1 . 

3. Global spectral fits with a plane-shock model imply 
strong overabundances of the elements we detect, con- 
sistent with their interpretation as ejecta dominated by 
intermediate-mass elements and Fe. 

4. We detect for the first time with high significance a line 
at 4.1 keV which we attribute to 44 Sc. The line strength 
of about 1.2 x 10~ 6 ph cm" 2 s" 1 implies a mass of Ti 
of 3.3(0.95,6.5) x 1O" 5 M , after correction for absorp- 
tion and assuming an age of 100 yr and a distance of 8.5 
kpc. 

5. We detect 44 Sc separately in the northern shell and the 
interior, with greater strength in the shell. 

We believe the discovery of 44 Sc emission to be the most 
important new result. A feature at 4. 1 ke V identified with 44 Sc 
has been reported previo u sly in G266.2-1.2 (RX J0 852.0- 
4622) by [T sunemi eTal] ([2000), Iyud in et al.| (pOTE) , and 
|Bamba et al.| ( |2005|>. flowever |Slane et al.| ( |2001| i and, more 
comprehensively, |Hiraga et~al . (2009) have failed to confirm 
these reports, at a level below the previous detection claims. It 
appears that Gl. 9+0.3 shows the first definite detection of this 
transition. Our line strength implies a flux in the 1.157 MeV 
line of about 1.2 x 10" 5 ph cm" 2 s" 1 , we ll below the COMP - 
TEL limit of about 2 x 10" 5 ph cm" 2 s" 1 ( |Dupraz et al.|l997| l. 



While our inferred 44 Ti mass is nominal ly in conflict with th e 
upper limit of 2 x 10" 5 M Q reported by Renaud et al. (2009), 
the errors are large enough to accommodate it. 

^Ti, a trace product of explosive nucleosynthesis, car- 
ries important inform ation ab out the det ails of the explosion 
(Diehl & Timmes 1998; The et al.|2006| l. It is primarily pro- 
duced in "a-rich freezeout" conditions, in which rapid cool- 
ing causes departures from nuclear statistical equilibrium and 
a high concentration of free a particles. In core-collapse 
(CC) SNe, its production is sensitive to the location of the 
mass cut (within which material becomes the neutron star), 
and in most models its production is correlated with that of 
56 Ni, and therefore with the supernova luminosity. Pr edicted 



yields r ange from 3 to 9 x 10 " M Q (see summary in The et 
al.|2006[ ). Traditional SN la models undergo a-rich freezeout 
under somewhat different conditions and do not produce large 
quantities of 44 Ti. However, recent simulations of off-center 
delayed-detonati on explosions have shown substantially in 
creased yields: 
of 1.6 x 10" 5 M ( 



Maeda et al. 



(2010 



0. 



compared to 2 x 



»)■ 
10 



quote a mass of ^Ti 
6 M for a centrally - 
6 M m for a central de- 



ignited pure deflagration and 3 x 10 iviq 
layed detonation. (For comparison, the benchmark W7 Type 
la model [Nomoto et al. 1984] give s 8 x 10~ 6 M o , while a se- 
ries of more recent simulations by Iwamoto et al. ( |1999 1 gives 
a range from (1 - 5) x 10" 5 M Q .) 

Galactic constraints on the supernova rate can be obtained 
from the observed abundance of 44 Ca, and the absence of ob- 
vious discrete sources of 1.157 MeV 7-ray emission, with 
the imp ortant exception of Cas A , the only firm detection 



to date (Iyudin et al. 1994 19991. The accompanying hard 



X-ray lines have also bee n seen from Cas A (Vink et al.| 



2001 



Renaud et al. l|2006| l, with an inferred mass of 44 Ti of 



about 2 x 1O" 4 M , more than expected f rom most CC models 
((6-14) x 1O~ 5 M ; [Timmes et al.|l996|, thoug h asymmetric 
models can do better (e.g., |Nagataki et al.|1998[ ). 

Since the 4. 1 keV inner-shell transition of 44 Sc formed from 
electron capture in 44 Ti does not require that the scandium be 
ionized, we are sensitive to both shocked and unshocked ma- 
terial. This is evident in the interior spectrum in Figure [3] in 
which the 44 Sc line is as prominent as in the north rim, but the 
other lines are weaker. The greater abundance of 44 Sc in the 
north, then, represents a true spatial distribution - surprising, 
considering that the 44 Sc is expected to be produced in the 
neighborhood of the Fe-peak elements, that is, the innermost 
ejecta. These properties suggest that the explosion was sub- 
stantially asymmetric, perhaps consistent with the greater pro- 
duction of T i found in asymmetric explosion models such as 
Maed a"etaLl ( |2010[ ). While He-shell denotations that trigger 
sub-Chandrasekhar thermonuclear explosions can also be in- 
voked to explain the presence of asymmetrically-distributed 
44 Ti and Fe in the outermost ejecta layers, too much 44 Ti is 
produced even in sub-Chandr asekhar ex plosion models with 
minimum He shell masses ( |Fink et al.|20 10). A substantially 
longer observation should allow better spatial localization of 
the 44 Sc emission, and an unprecedented view into supernova 
nucleosynthesis. 

Our discovery of these spectral lines, and the 44 Sc line 
in particular, remains consistent with the possibility that 
Gl.9+0.3 resulted from a Type la event, a possibility we find 
increasingly likely. The bilaterally symmetric synchrotron- 
X-ray morphology, the extremely high shock velocities we 
both infer and measure from line broadening, and the promi- 
nent Fe K emission, all support a Type la origin. We have not 
yet ruled out a CC origin, however; while strong iron emission 
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Table 1 

Spectral Fits 



Region Nh 

(10 22 cm- 2 ) 



F5-10 keV 

(1CT 13 ergs crrT 2 s" 1 ) 



Vshell 



Line energies (keV) and line strengths (10 7 ph cm 2 s ') 



(10 4 kms-') 


Silicon 


Sulfur 


Argon 


Calcium 


Scandium 


Iron 


1.4 


1.1 


1.83 


2.35 


2.99 


3.72 


4.08 


6.49 


(0.9, 1.9) 


< 1.9 


(1.78, 1.88) 


(2.31,2.40) 


(2.92, 3.04) 


(3.59, 3.87) 


(4.00, 4.15) 


(6.38,6.61) 






4.4 


10.4 


7.7 


5.2 


10 


18 






(2.1,7.0) 


(5.9, 15.6) 


(3.1, 13) 


(1.2, 11) 


(4.0, 19) 


(11, 27) 


2.1 


0.9 


1.83 


2.32 


2.96 


3.74 


4.09 


6.4 


(1.2,3.0) 


< 1.6 


(1.72, 1.95) 


(2.17,2.46) 


(2.87, 3.05) 


(3.52, 3.94) 


(4.02, 4.17) 


(6.0, 6.8) 






1.2 


3.1 


5.9 


1.9 


4.6 


3.0 






(0.1,4.0) 


(0.4, 8.3) 


(1.5, 13) 


(0.2, 6.0) 


(1.3,9.5) 


(0.4, 8.4) 



North rim 6.88 2.40 3.00 

(6.72,7.05) (2.29,2.51) (2.80,3.22) 



Interior 6.89 2.71 1.35 

(6.72,7.06) (2.54,2.86) (1.22,1.51) 



Note. — Line energies (strengths) are in rows 1 and 5 (3 and 7), with 90% confidence limits listed in adjacent rows. 



is unusual in CC remnants, Cas A is a counterexample, and it 
is possible that a CC origin can explain the other features as 
well. What is certainly clear is that G 1.9+0. 3 still has a great 
deal to teach us about the evolution of very young SNRs and 
about the supernovae that produced them. 

This work was supported by NASA through Chandra Gen- 
eral Observer Program grant SAO GO6-7059X. 
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